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ABSTRACT 

We explore the utility of future photometric redshift imaging surveys for delineating 
the large-scale structure of the Universe, and assess the resulting constraints on the 
cosmological model. We perform two complementary types of analysis: 

(1) We quantify the statistical confidence and the accuracy with which such sur- 
veys will be able to detect and measure characteristic features in the clustering power 
spectrum such as the acoustic oscillations and the turnover, in a 'model-independent' 
fashion. We show for example that a 10,000 deg 2 imaging survey with depth r = 22.5 
and photometric redshift accuracy Sz/(1 + z) — 0.03 will detect the acoustic oscilla- 
tions with 99.9 per cent confidence, measuring the associated preferred cosmological 
scale with 2 per cent precision. Such a survey will also detect the turnover with 95 per 
cent confidence, determining the corresponding scale with 20 per cent accuracy. 

(2) By assuming a ACDM model power spectrum we calculate the confidence with 
which a non-zero baryon fraction can be deduced from such future galaxy surveys. We 
quantify 'wiggle detection' by calculating the number of standard deviations by which 
the baryon fraction is measured, after marginalizing over the shape parameter. This is 
typically a factor of four more significant (in terms of number of standard deviations) 
than the above 'model-independent' result. 

For both analyses we quantify the variation of the results with magnitude depth 
and photometric redshift precision, and discuss the prospects for obtaining the required 
performance with realistic future surveys. We conclude that the precision with which 
the clustering pattern may be inferred from future photometric redshift surveys will be 
competitive with contemporaneous spectroscopic redshift surveys, assuming that sys- 
tematic effects can be controlled. We find that for equivalent wiggle-detection power, a 
photometric redshift survey requires an area approximately 12(5 z/(l + z)/0.03) times 
larger than a spectroscopic survey, for a given magnitude limit. We also note that 
an analysis of Luminous Red Galaxies in the Sloan Digital Sky Survey may yield a 
marginal detection of acoustic oscillations in the imaging survey, in addition to that 
recently reported for the spectroscopic component. 

Key words: large-scale structure of Universe - cosmological parameters - surveys 



1 INTRODUCTION 



be discriminated by accurate measurements of the shape of 
the primordial power spectrum of mass fluctuations? 



Today's most pressing cosmological questions demand the 
construction of galaxy surveys of unprecedented depth and 
volume. Such questions include: is the accelerating rate of 
cosmic expansion driven by Einstein's cosmological constant 
or some different form of 'dark energy'? What are the prop- 
erties of this dark energy? Can competing models of inflation 



Galaxy surveys delineate the large-scale structure of 
the Universe and thereby provide a powerful and indepen- 
dent constraint on the cosmological model. The currently- 
favoured 'concordance model' - in which w 70 per cent of the 
energy density of today's Universe is resident in a relatively 
unclustered form known as 'dark energy' - is evidenced by a 
combination of observations of the Cosmic Microwave Back- 
ground (e.g. Spergel et al. 2003) with either those of galaxy 
clustering (e.g. Percival et al. 2001) or of high-redshift super- 
novae (e.g. Riess et al. 1998; Perlmutter et al. 1999). Either 
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pair of these independent datasets are required to break the 
degeneracies between model parameters and render a unique 
cosmology. 

According to standard cosmological theory, if the linear- 
regime clustering power spectrum is measured with suffi- 
cient precision then it will no longer appear smooth and 
monotonic: specific features and modulations will become 
apparent. Two such attributes are predicted: firstly, a series 
of acoustic oscillations - sinusoidal modulations in power 
as a function of scale imprinted in the baryonic component 
before recombination (Peebles & Yu f970; Hu & Sugiyama 
f996) - and secondly, a turnover - a broad maximum in clus- 
tering power on large scales originating from the radiation- 
dominated epoch. These features encode characteristic cos- 
mological scales that can be extracted from the observa- 
tions, greatly improving constraints upon cosmological mod- 
els (e.g. Blake & Glazebrook 2003; Seo & Eisenstein 2003). 
Moreover, other currently-unknown modulations in power 
(e.g. signatures of inflation) may be discovered when the 
clustering pattern is examined with sufficiently high preci- 
sion (e.g. Martin & Ringeval 2004). 

Very recently, the acoustic signature has been convinc- 
ingly identified for the first time in the clustering pattern 
of Luminous Red Galaxies in the Sloan Digital Sky Sur- 
vey (SDSS; Eisenstein et al. 2005). The 2dF Galaxy Red- 
shift Survey has produced consistent measurements (Cole 
et al. 2005). These results confirm previous tantalizing hints 
(e.g. Miller, Nichol & Chen 2002; Percival et al. 2001). The 
challenge now is to make more accurate measurements at 
different redshifts, using these features to further constrain 
the cosmological parameters, in particular the dark energy 
model. At low redshift the available volume is limited: the 
effect of cosmic variance is significant. Therefore such sur- 
veys are insensitive to clustering modes on very large scales 
and are hampered by non-linear growth of structure on small 
scales. Higher-redshift large-scale surveys are consequently 
required to map greater cosmic volumes: tracing clustering 
modes with longer wavelengths and additionally unveiling 
the pattern of linear clustering to significantly smaller scales. 

The high-redshift spectroscopic surveys currently being 
executed (e.g. DEEP2, Davis et al. 2003; VVDS, Le Fevre 
et al. 2003) cover solid angles of ~ 10 deg 2 , which are insuf- 
ficient for detecting the predicted features in the clustering 
power spectrum. Such projects are fundamentally limited 
by existing instrumentation, being performed by spectro- 
graphs with relatively small fields-of-view (^ 10 — 20') and 
restricted (albeit impressive) multi-object capabilities. Some 
proposed new instrumentation addresses this difficulty (e.g. 
the KAOS project, Barden et al. 2004), permitting spectro- 
scopic exposures over single fields of ^ 1 deg 2 using £s 5000 
fibres. However, these projects will take many years to reach 
completion. 

In this paper we consider the role that photometric 
redshift catalogues derived from deep imaging surveys could 
play in addressing the scientific goals outlined above. Ex- 
tensive imaging surveys (covering ~ 10,000 deg 2 ) to rea- 
sonable depths (r ~ 22) are ongoing (e.g. SDSS); the im- 
plied redshift distributions map galaxies over cosmic dis- 
tances to z m 1 with sufficient number density that clus- 
tering measurements are limited by cosmic variance rather 
than by shot noise. Future deeper imaging surveys (e.g. Pan- 
STARRS, Kaiser et al. 2000; CTIO Dark Energy Survey, 



http://cosmology.astro.uiuc.edu/DES; LSST, Tyson et 
al. 2002) are being planned to address a host of scientific 
questions including in particular weak gravitational lensing. 
We argue that such surveys will also provide powerful mea- 
surements of features in the galaxy clustering pattern. 

The utility of photometric redshifts - derived from 
broadband galaxy colours rather than from spectra - has 
been well-established, with many different techniques being 
successfully utilized. The simplest method involves the fit- 
ting of model spectral templates (e.g. Bolzonella, Miralles 
& Pcllo 2000). Other approaches use spectroscopic 'training 
sets' to calibrate the photometric redshifts via an empiri- 
cal polynomial of colour terms (Connolly et al. 1995) or an 
artificial neural network (Firth, Lahav & Somerville 2003). 
The precision Sz with which galaxy redshifts (and therefore 
radial distances) may be determined varies with the method 
and filter set used, together with the galaxy type, magnitude 
and redshift, but at best is currently <to = Sz/(1 + z) ~ 0.03 
(e.g. COMBO-17, Wolf et al. 2003). For the SDSS imaging 
component, the rms photometric redshift accuracy varies 
from Sz ~ 0.03 for bright galaxies with r < 18 to Sz ~ 0.1 
for magnitudes r ~ 21 (Csabai et al. 2003). 

The blurring of large-scale structure in the radial direc- 
tion due to the photometric redshift error degrades measure- 
ments of the clustering pattern. However, on physical scales 
larger than that implied by the redshift error, the informa- 
tion is preserved. Moreover, on smaller scales the tangen- 
tial information always survives, and the vast area which 
may be readily covered by an imaging survey can poten- 
tially provide more independent structure modes on a given 
scale than those yielded by a fully spectroscopic survey of a 
smaller solid angle, implying very competitive cosmological 
constraints. Photometric redshifts have already been used 
to construct volume-limited samples of low-redshift galax- 
ies and measure their angular clustering properties (Bu- 
davari et al. 2003; see also Meiksin, White & Peacock 1999, 
Cooray et al. 2001). The cosmological parameter constraints 
resulting from future photometric redshift imaging surveys 
have been simulated by Seo & Eisenstein (2003) ; Amendola, 
Quercellini & Giallongo (2004) and Dolney, Jain & Takada 
(2004). 

In this study we use a Monte Carlo approach to model 
the galaxy power spectra resulting from a host of simulated 
photometric redshift surveys as a function of the limiting 
magnitude of the initial imaging and the accuracy of the 
derived photometric redshift. Our simulation methodology 
is described in Section 2, where first results for the accu- 
racy of power spectrum measurements are presented. We 
infer constraints on the cosmological model using two com- 
plementary methods with very different prior assumptions. 
Firstly, in Sections 3 and 4 we discuss in detail the re- 
sulting confidence of detection of the acoustic oscillations 
and power spectrum turnover, respectively, and the accu- 
racy with which the associated characteristic cosmological 
scales may be extracted. In these analyses we make minimal 
assumptions, purely concerning ourselves with the statistical 
detection of power spectrum features relative to a smooth 
monotonic fit. Secondly, in Section 5 we use the full power 
spectrum shape information in conjunction with theoretical 
fitting formulae to compute constraints on the basic param- 
eters of the cosmological model, in particular the baryon 
fraction f2b/n m and the running of the spectral index of the 
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primordial power spectrum n run - In all cases we compare our 
results to those deduced from spectroscopic redshift surveys. 
We evaluate the effect of the approximations of our method- 
ology in Section 6, in particular considering a wider range of 
photometric redshift error distributions. Finally, in Section 
7 we outline the prospects for obtaining the requisite imag- 
ing depth and photometric redshift accuracy using realistic 
future surveys. 



2 MODELLING THE POWER SPECTRUM OF 
PHOTOMETRIC REDSHIFT SURVEYS 

2.1 Method summary 

Our methodology for simulating the large-scale structure of 
future galaxy surveys is to generate many 'Monte Carlo' re- 
alizations of galaxy distributions from a fiducial power spec- 
trum. An observed power spectrum is measured for each re- 
alization separately, using techniques similar to those which 
would be employed for real survey data. The resulting en- 
semble of observed power spectra can then be used to quan- 
tify the error distribution in derived quantities, without any 
need to approximate the likelihood surface by techniques 
such as Fisher matrices (e.g. Seo & Eisenstein 2003; Amen- 
dola et al. 2004; Dolney et al. 2004). For example, the stan- 
dard deviation in the measurement of the power spectrum 
P(k) in a given bin around scale k follows from the scat- 
ter in the recovered values of P (k) in that bin across the 
realizations, without the need for analytic approximations. 

Our procedure for modelling photometric redshifts is 
to convolve the simulated galaxy distributions in the ra- 
dial direction with a photometric redshift error function (in 
the fiducial case, a Gaussian with width a x in real space). 
The measured power spectrum is derived by computing 
the Fourier transform of the whole survey box, then dis- 
carding small-scale radial Fourier modes with wavenumbers 
fcrad ^ (which contribute no signal due to the radial 

smearing) . Note that the resulting number of useful Fourier 
structure modes is identical to that obtained if the survey 
box is instead split into many independent slabs of width a x 
and a purely angular power spectrum is measured for each 
slice. 

In particular, we wish to assess the confidence with 
which we can detect specific features in the clustering power 
spectrum such as the acoustic oscillations and the 'turnover'. 
These features can be modelled by simple empirical formulae 
which can be fitted to each measured power spectrum real- 
ization, and the best-fit \ 2 statistic calculated. The resulting 
best-fit x 2 can be compared with that of a smooth (feature- 
less) power spectrum fit, resulting in a relative probability 
of feature detection for each realization. The distribution of 
the relative probabilities across the Monte Carlo realizations 
enables a very realistic assessment of the efficacy of future 
surveys across the statistical ensemble of possible universes. 

Furthermore, we are interested in recovering character- 
istic scales from these features in the power spectrum. These 
scales can be encoded into our empirical fitting formulae; the 
distribution of best-fitting values of these scales across the 
realizations is indicative of the realistic accuracy with which 
it is possible to measure them with the simulated survey. 

In order to perform our simulations we must also adopt 



a fiducial set of cosmological parameters, which determine 
both the cosmic volume mapped by a given survey and the 
fiducial power spectrum (via the fitting formulae of Eisen- 
stein & Hu 1998). 

We characterize a photometric redshift imaging survey 
using two parameters: 

• The photometric redshift error distribution (in the sim- 
plest case, a Gaussian distribution with standard deviation 
Sz), which controls the 'smearing' of the underlying large- 
scale structure in the radial direction. 

• The threshold apparent magnitude of the input imaging 
catalogue, mi im , which determines the redshift distribution 
dN/dz of the 'unsmeared' galaxy distribution, i.e. the radial 
depth of the survey. This magnitude limit is defined using 
the SDSS r-band filter. 

In this paper we present the results of simulations of a range 
of photometric redshift surveys as a function of these two 
parameters. We assume the survey area in all cases is 10,000 
deg 2 . 

In the following Section we provide a detailed account 
of the assumptions and method we used to simulate the 
observed power spectra. For analyses of future spectroscopic 
redshift surveys using a similar method we refer the reader to 
Blake & Glazebrook (2003) and Glazebrook & Blake (2005). 



2.2 Detailed simulation methodology 

(i) A fiducial cosmology is chosen for the simulation. Un- 
less otherwise stated we assumed a flat ACDM Universe with 
parameters fi m = 0.3, J1a = 0.7, h — 0.7 and fib/f2 m = 0.15. 

(ii) A limiting apparent magnitude mn m is assumed 
for the imaging survey. A model redshift distribution 
dN/dz{m\\- m ) is determined, as described in Section 2.3. 

(iii) A survey redshift range (z m i n ,z max ) and solid angle 
An is specified. We assumed An — 10,000 deg 2 for our imag- 
ing surveys. The chosen redshift interval depends on mn m 
as described in Section 2.3 (and listed in Table 1). 

(iv) We performed our simulations using a 'flat-sky ap- 
proximation' for computational convenience (this approxi- 
mation has a negligible effect on our results as discussed 
in Section 6). A cuboid with sides of co-moving lengths 
(L x ,L y ,L z ) is created, possessing a volume equal to that 
enclosed by the survey cone. We take the x-axis as the ra- 
dial direction. The length L x is the co-moving distance be- 
tween redshifts z m i n and z max , and the other dimensions are 
determined by stipulating L y = L z (although the results 
are independent of the ratio L y /L z , assuming that both of 
these dimensions are large enough to imply a sensitivity to 
structural modes with scales contributing to the acoustic 
oscillations) . 

(v) A model linear theory matter power spectrum 
fmass(fc,z = 0) is computed for the chosen parameters 
(f2 m ,fib,/i) from the fitting formula of Eisenstein & Hu 
(1998), assuming a z — normalization cr 8 = 1 and a pri- 
mordial power-law slope n = 1. The survey is assumed to 
have an 'effective' redshift z e g — (z m in + Zmax)/2. The power 
spectrum is scaled to this redshift using a standard ACDM 
growth factor: 

Pgal(fc, Zcff) = -Pmass(fc, 0) D 1 (z c f{) 2 fo 2 (1) 
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where we use the Carroll, Press & Turner (1992) approx- 
imation for Di(z) and a constant linear bias factor b for 
the clustering of galaxies with respect to matter. The value 
b = 1 is assumed for our surveys, unless otherwise stated. 
We neglect any evolution of clustering across the depth of 
the survey box. This approximation is discussed in Section 6; 
we note that the clustering amplitude of galaxies is known 
to evolve much less rapidly with redshift than that of the 
underlying mass fluctuations (indicating an evolution of the 
galaxy bias parameter in the opposite sense). 

(vi) The location in fc-space of the transition between the 
linear and non- linear regimes of gravitational clustering, kn n , 
is determined from P ga i(fc) in a conservative manner (see 
Blake & Glazebrook 2003, Figure 1). We only measure power 
spectra in the linear regime, i.e. for scales k < kn n . 

(vii) A set of Monte Carlo realizations (numbering 400 for 
all simulations presented here) is then performed to generate 
many different galaxy distributions consistent with P ga i(fc), 
as described in steps (viii) and (ix). 

(viii) A cuboid of Fourier coefficients is constructed with 
grid lines set by dki = 2ir/Li, with a Gaussian distribution of 
amplitudes determined from P ga i(fc), and with randomized 
phases. The gridding is sufficiently fine that the Nyquist 
frequencies in all directions are significantly greater than 

fclin- 

(ix) The Fourier cuboid is FFTed to determine the den- 
sity field in the real-space box. The result is modulated by 
the survey redshift distribution dN/dz. This observed den- 
sity field is then Poisson sampled to determine the number 
of galaxies in each grid cell. 

(x) A photometric redshift error distribution is assumed. 
For our main set of simulations, we modelled this function 
as a Gaussian distribution, such that the radial co-moving 
co-ordinate x of each galaxy was smeared by an amount 8x 
sampled from a probability distribution 



f(5x) oc exp 



2 \a x ) 



(2) 



In practice we specified a redshift error parameter a and 
derived a x in accordance with the equation 



dx c 
ax=5z _ {z = Zes) = ao{1 + Zes) __ 



(3) 



where H(z) is the value of the Hubble constant measured by 
an observer at redshift z. Equation 3 encodes the expected 
zeroth-order dependence of photometric redshift precision, 
5 z oc (1 + z), originating from the stretching of galaxy spec- 
tra with redshift for a filter system with constant spectral 
resolution AA/A. We assess the effect of more complex pho- 
tometric redshift error distributions than Equation 2 in Sec- 
tion 6.1. 

(xi) The galaxy number distribution is 'smeared' along 
the x— (radial) direction in accordance with the photomet- 
ric redshift error function, taking pixelization effects into 
account. The resulting distribution is our simulated photo- 
metric redshift survey. We note that our simple photometric- 
redshift error model represents a convolution of the 'un- 
smeared' galaxy number distribution with the error function 
fix) (Equation 2). In this case, according to the convolu- 
tion theorem, the resulting power spectrum signal is damped 
along the radial direction: 

P(k x ,k y ,k z ) -> P(k x ,k y ,k z ) X exp [-(k x a x ) 2 ] (4) 



where g(k x ) — exp [— (k x a x ) 2 ] is the square of the Fourier 
transform of f(x). 

(xii) The power spectrum of the resulting distribution is 
measured using standard estimation tools: essentially this 
involves taking the Fourier transform of the density field, 
subtracting that of the survey window function, and bin- 
ning up the resulting modes in fc-space (see e.g. Hoyle et al. 
2002; note that we do not use the optimal-weighting tech- 
nique presented by Feldman, Kaiser & Peacock (1994) be- 
cause this does not represent a simple convolution of the 
density field and consequently step (xiii) below would not 
be possible). Power spectrum modes in Fourier space are 
divided into bins of total k, up to a maximum of fcij n . We 
only include modes with a value of | k x | less than a maximum 
This is because in accordance with Equation 4, the 
photometric-redshift smearing damps the clustering signal 
along the radial direction such that modes with high values 
of Ifc^l contribute only noise. Hence the dominant contribu- 
tion to power spectrum bins with k > fe 1)Illal originates from 
tangential Fourier modes with k x ~ and yjky + k'i ~ k. 
The value of fc^max is determined by the equation 



2/a x 



(5) 



where the coefficient of 2 was determined by experiment to 
be optimal for the surveys presented here. Use instead of a 
coefficient of 1.5 does not change the results significantly, 
but 1.0 is sub-optimal. 

(xiii) The measured power spectrum P(k) is 'undamped' 
by dividing by a function /damp {k). This 'damping function' 
was determined by binning the expression exp [— (k x a x ) 2 ] 
(from Equation 4) as a function of total k as described in 
step (xii). 

(xiv) An error bar is assigned to each power spectrum bin 
using the variance measured over the Monte Carlo realiza- 
tions. 

We do not incorporate redshift- space distortions into our 
simulations because the implied radial smearing due to pecu- 
liar velocities is much less than that due to the photometric 
redshift error. 



2.3 Modelling the redshift distribution 

In order to model the 'unsmeared' survey redshift distri- 
bution as a function of the limiting apparent magnitude of 
the imaging survey mii m , we used the luminosity functions 
derived from the COMBO-17 survey (Wolf et al. 2003). Ta- 
ble A. 2 of Wolf et al. (2003) lists Schechter function pa- 
rameters in redshift slices of width Az — 0.2 in the range 
0.2 < 2 < 1.2 for the SDSS r filter. In the regime z < 0.2 
we applied the locally-determined SDSS luminosity function 
(Blanton et al. 2003). For a given threshold apparent mag- 
nitude rii m , we converted these luminosity functions into a 
redshift distribution dN/dz for each redshift slice, fitting the 
overall result with a simple model parameterized by a char- 
acteristic redshift zo and an overall surface density Eo (in 
deg- 2 ): 



dN 

dz 



3z 2 

^ = Eo x 2F exp 



(if 



(6) 
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Figure 1. Fractional error in P(k) against k for examples of photometric redshift imaging surveys. The left-hand panel illustrates the 
variation of the fractional power spectrum error as a function of r\i m for <ro = 0.03; the right-hand panel displays the dependence on 
fTQ for ri im = 22. We bin the power spectra in intervals of Ak = 0.005 h Mpc -1 . Simulated results for the SDSS main spectroscopic 
survey (dashed line) and SDSS LRG spectroscopic survey (dotted line) are overplotted (see Section 2.4 for details of our models of these 
surveys). Considering the photometric redshift surveys, the fractional power spectrum precision improves with increasing survey depth 
(left-hand panel) owing to the greater cosmic volume mapped and the consequently higher density-of-states in fc-space. The performance 
degrades with increasing <ro (right-hand panel) owing to the decreasing width of the 'undamped' slab in fc-space (i.e. the decreasing value 
of &a;,max), although on very large scales {k < k x ,max) all Fourier modes are retained and performance is unaffected. 



(e.g. Baugh & Efstathiou 1993). The values of the fitted 
parameters are displayed in Table 1. We used model K- 
corrections averaged over different galaxy spectral types. 

For each apparent magnitude limit we selected a red- 
shift interval (z m i n , z ma i) for the simulation. For all but the 
shallowest surveys we set z m i n = 0.2, the results are insen- 
sitive to this choice because there is minimal volume con- 
tained by lower redshifts. As the value of z max increases, the 
variance in the recovered power spectrum is determined by 
a balance between two conflicting effects: increasing survey 
volume (i.e. decreasing cosmic variance) and decreasing av- 
erage number density (i.e. increasing shot noise) owing to 
the fixed magnitude threshold. We determined the optimal 
value of Zmax for each magnitude threshold by experimenting 
to determine the most accurate measurement of the acoustic 
oscillations (see Section 3). Our chosen ranges are listed in 
Table 1. 

We note that the optimal value of z max for measuring 
power spectrum modes around the turnover is marginally 
higher than that for detecting the acoustic oscillations, be- 
cause in the former case the power spectrum amplitude is 
at a maximum, implying a lower required galaxy number 
density for suppressing shot noise. We always use the more 
conservative values of z max in Table 1, but this does not 
change our results significantly. 



Table 1. Input parameters for simulated galaxy redshift surveys 
as a function of limiting apparent magnitude in the SDSS r filter, 
Him- The 'unsmeared' redshift distribution is specified by the val- 
ues of 20 and So in accordance with Equation 6. The minimum 
and maximum redshifts of the simulated survey, z m i n and z ma x, 
are also listed. 



Him 


zo 


So (deg- 2 ) 


•2m in 


^max 


18.0 


0.1 


120 


0.1 


0.4 


18.5 


0.12 


230 


0.1 


0.5 


19.0 


0.14 


410 


0.1 


0.6 


19.5 


0.16 


710 


0.1 


0.6 


20.0 


0.18 


1200 


0.2 


0.7 


20.5 


0.2 


2000 


0.2 


0.7 


21.0 


0.22 


3200 


0.2 


0.8 


21.5 


0.24 


4900 


0.2 


0.8 


22.0 


0.27 


7500 


0.2 


0.9 


22.5 


0.3 


11100 


0.2 


1.0 


23.0 


0.33 


16300 


0.2 


1.2 


23.5 


0.36 


24000 


0.2 


1.3 


24.0 


0.39 


35400 


0.2 


1.4 



We note that all magnitudes referred to here are total 
galaxy magnitudes, and we neglect all incompleteness effects 
associated with surface brightness. These are expected to be 
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minimal at high redshift owing to the decreasing apparent 
size of galaxies. We also neglect the fact that imaging surveys 
to depths rii m > 23 may be used to select additional galaxy 
populations at z ^ 3 using the Lyman Break technique. 

2.4 Fractional errors in the power spectrum 

We simulated photometric redshift surveys with limiting 
magnitudes varying from rn m = 20 to 24 and photomet- 
ric redshift error parameters in the interval cro = 0.01 to 

0. 05. These ranges were chosen to encompass the current 
state-of-the-art and realistic future improvements. Achiev- 
ing precision cro = 0.01 would require observations with 
many narrow-band filters, probably encompassing the near 
infra-red. 

As a first application of our method, Figure 1 displays 
the fractional standard deviation in the power spectrum 
measurement as a function of scale (averaged over the Monte 
Carlo realizations) for some example photometric redshift 
surveys. The measured power spectrum modes are averaged 
over angles in bins of width Afc = 0.005 h Mpc -1 . 

The scaling of the resulting power spectrum errors as a 
function of fc can be understood simply by counting the num- 
ber of Fourier modes m within each bin (the errors SP scal- 
ing as 1 /^Jrn). For a bin with fc 3> fcc.max, these modes are lo- 
cated approximately within an cylindrical annulus in Fourier 
space of radius k, thickness Afc and depth fc x , ma x. This 
amounts to a volume in fc-space equal to 27rfc x k x .max x Afc, 

1. e. SP oc fc -1 / 2 . This contrasts with a fully spectroscopic 
survey, for which the relevant Fourier modes for a scale fc 
reside within a spherical shell, such that SP oc fc -1 . Al- 
though a photometric redshift survey maps out a reduced 
volume of Fourier space, the larger density-of-states (owing 
to the increased cosmic volume probed) can still result in a 
more accurate measurement of the galaxy power spectrum. 

In order to illustrate this point, we compare the power 
spectrum accuracies for our simulated photometric redshift 
surveys with those expected for the SDSS spectroscopic sur- 
veys (both the main survey and the Luminous Red Galaxy 
(LRG) survey). We created Monte Carlo power spectrum 
realizations for these SDSS spectroscopic surveys using the 
methodology of Blake & Glazebrook (2003), which is very 
similar to that presented in Section 2 above. The main differ- 
ences are that no photometric redshift smearing is applied, 
and thus a conical geometry may be employed rather than a 
flat-sky approximation (although as discussed in Section 6, 
this makes a negligible difference to the results). We mod- 
elled the SDSS main spectroscopic survey using the redshift 
distribution 



dN 2 
— — oc z exp 

dz 



(0.O55) 



(7) 



with a total surface density equal to Eo = 70.7 deg -2 . This 
model constitutes a good fit to the relevant luminosity func- 
tion (Blanton et al. 2003). The redshift interval of the simu- 
lation was < z < 0.25. We approximated the LRG spectro- 
scopic survey redshift distribution using the Gaussian func- 
tion 



dN 

dz 



oc exp 



/ z- 0.375 \ 
V 0.065 ) 



(8) 



and a total surface density Eo = 17.2 deg 2 , which provides 



a reasonable fit to the radial selection function discussed by 
Eisenstein et al. (2001). The redshift interval of the simu- 
lation was 0.3 < z < 0.45 and LRGs are assigned a linear 
bias factor 6 = 2 (see Equation 1). For both spectroscopic 
surveys we assumed an areal coverage of 10,000 deg 2 , the 
same as for the simulated photometric redshift surveys. 

As illustrated by Figure 1, in the turnover regime 
(fc < 0.02 h Mpc -1 ) the photometric redshift surveys always 
yield more large-scale modes than the SDSS spectroscopic 
surveys owing to the larger cosmic volume mapped and the 
fact that the wavelengths of these modes significantly ex- 
ceed the length-scale of photometric-redshift radial smear- 
ing. In the acoustic oscillations regime (fc > 0.05 ft Mpc -1 ), 
a 10,000 deg photometric redshift survey out-performs the 
SDSS LRG spectroscopic survey provided that 00 < 0.03 and 
riim > 22.5. 



3 MEASURING THE ACOUSTIC 

OSCILLATIONS USING PHOTOMETRIC 
REDSHIFT SURVEYS 

The clustering power spectrum on intermediate scales 
(0.05 < fc < 0.3/i Mpc -1 ) contains small-amplitude (~ 5 
per cent) modulations known as 'acoustic oscillations' (Pee- 
bles & Yu 1970; Hu & Sugiyama 1996; Eisenstein & Hu 

1998) . These fluctuations in power have an identical physi- 
cal origin to those observed in the CMB, namely, oscillations 
in the photon-baryon fluid before recombination. 

There has been considerable recent interest in exploit- 
ing these acoustic features as an accurate and clean probe 
of the cosmological model (Blake & Glazebrook 2003; Seo & 
Eisenstein 2003; Linder 2003; Hu & Haiman 2003). The ap- 
proximately sinusoidal fluctuations in power encode a char- 
acteristic scale - the sound horizon at recombination - which 
can be measured from the CMB. This scale can then act as 
a standard cosmological ruler (Eisenstein, Hu & Tegmark 

1999) : its recovered value from a galaxy redshift survey de- 
pends on the assumed cosmological parameters, in partic- 
ular the dark energy model, and may be used to constrain 
those parameters in a manner that is probably significantly 
less sensitive to systematic error than other probes (Blake 
& Glazebrook 2003). 

Very recently, analysis of the clustering pattern of SDSS 
Luminous Red Galaxies at z « 0.35 has yielded the first 
convincing detection of the acoustic signal and application 
of the standard ruler (Eisenstein et al. 2005) . Although this 
survey does not have sufficient redshift reach to constrain 
dark energy models, this result is an important validation 
of the technique. Indeed, detection of these acoustic features 
represents a fundamental test of the paradigm of the origin 
of galaxies in the fluctuations observed in the CMB. 

Utilization of the acoustic oscillations to measure the 
properties of dark energy demands new galaxy surveys of un- 
precedented depth and volume (Blake & Glazebrook 2003; 
Seo & Eisenstein 2003). Given the current availability of 
large-scale imaging surveys such as the SDSS, and the an- 
ticipated wait of several years for commencement of projects 
with sufficiently capable spectroscopic facilities able to sur- 
vey ~ 10 6 objects over ~ 1000 deg 2 (such as the KAOS 
proposal), it is timely to evaluate the role photometric red- 
shift surveys could play in the detection and measurement 
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Figure 2. Mean values of the power spectrum divided by the smooth reference spectrum, P(k)/P Ia {, averaged over the Monte Carlo 
realizations for various photometric redshift surveys, together with the standard deviation in each bin. The left-hand panel illustrates 
the variation of the observed P(k)/P tB f as a function of r Um for ctq = 0.03; the right-hand panel displays the dependence on <tq for 
ry lm = 22. Results for different surveys are offset for clarity. The power spectra are plotted for wavelengths larger than the linear/non- 
linear transition scale ku n . This scale varies with the threshold magnitude rn m because it is evaluated at the effective redshift of the 
survey z e ff = (z m i n + z max )/2 (see Table 1). In all cases we use a survey area of 10,000 deg 2 . 



of the acoustic oscillations. Furthermore, recently-developed 
novel photometric-redshift techniques such as those utilizing 
artificial neural networks should prove extremely useful in 
this regard (e.g. Collister & Lahav 2004). 

The constraints on the cosmological model yielded by 
acoustic oscillations in future photometric redshift surveys 
have been discussed by Seo & Eisenstein (2003); Amendola 
et al. (2004) and Dolney et al. (2004). Here we take a dif- 
ferent but complementary approach. Firstly, this previous 
work deduced cosmological parameter constraints using a 
Fisher matrix approach which provides the minimum pos- 
sible errors for an unbiased estimate of a given parameter, 
based upon the curvature of the likelihood surface near the 
fiducial model. In the present study we instead use Monte 
Carlo techniques, which make a closer connection with the 
analysis methods that would be used for real data and can 
probe more realistic non-parabolic likelihood surfaces. Sec- 
ondly, we give detailed consideration to the statistical confi- 
dence of detection of the relevant power spectrum features, 
carefully separating this information from that contained in 
the overall shape of the power spectrum, which may be sub- 
ject to additional systematic distortions, as discussed below. 
Thirdly, by treating a wide grid of potential photometric 
redshift surveys varying both the redshift accuracy and the 
limiting magnitude, we can make a direct connection with 
the performance of current and future experiments. 

The comparison of photometric redshift and spectro- 



scopic redshift surveys has already been discussed in detail 
by Blake & Glazebrook (2003) and Seo & Eisenstein (2003). 
To summarize the relevant points of these two papers: 

• As discussed in Section 2.2, Fourier modes with values 
of k x ^> 1 ja x (where x is the radial axis) contribute noise. A 
photometric redshift survey therefore requires significantly 
more sky area than a spectroscopic redshift survey of similar 
depth to yield the same number of Fourier modes in a given 
power spectrum bin with scale k ^> \ ja x . 

• Fourier modes with usable signal-to-noise ratios are 
largely tangential (k x w 0). Consequently, in the case of 
a photometric redshift survey, we are only able to apply 
the standard ruler represented by the acoustic oscillations 
in the tangential direction, constraining the co-ordinate dis- 
tance x(z) to the effective redshift of the survey. We lose 
the capacity of a spectroscopic redshift survey to apply the 
ruler radially, measuring dx/dz (or equivalently the Hub- 
ble constant at redshift z), which yields powerful additional 
constraints on the dark energy model. 

In this Section we present a series of simulations ad- 
dressing the issues of the confidence and accuracy of de- 
tection of acoustic oscillations as a function of photometric 
redshift error <to (as defined by Equation 2) and limiting 
apparent magnitude ra m of the imaging survey. We defer 
the questions of whether and how these requirements can 
be realized in realistic surveys to Section 7. We proceed in 
a model-independent fashion, quantifying the statistical sig- 
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nificance with which we can assert deviations from a feature- 
less monotonic function, without using any of the informa- 
tion contained in the power spectrum shape. In Section 5 we 
use a more model-dependent approach, combining the full 
power spectrum shape function with recent measurements 
of the CMB anisotropics to derive predicted constraints on 
the cosmological parameters. 

Our set of Monte-Carlo realized power spectra enables 
us to evaluate statistical questions of confidence and accu- 
racy over a realistic ensemble of Universes, without needing 
to approximate the statistical distributions or likelihood sur- 
faces, except that when converting values of the x 2 statistic 
to relative probabilities we implicitly assume that the errors 
in the measured power spectra are Gaussian, which agrees 
well with histograms obtained from the Monte Carlo real- 
izations. 



3.1 Confidence of detection of acoustic oscillations 

We note that confidence of detection can be defined in sev- 
eral different ways and depends strongly on statistical pri- 
ors. One approach to the data analysis would be to fit full 
ACDM transfer functions (e.g. the formulae of Eisenstein & 
Hu 1998) to the measured power spectra and thereby de- 
termine that baryonic models (containing acoustic features) 
provided a significantly better fit to the data than models 
with fib = 0. We argue in the current Section that this only 
partially constitutes a detection of acoustic oscillations, be- 
cause information contained in the shape of the power spec- 
trum is also constraining this fit. 

We adopted a conservative approach in which, prior 
to measuring the preferred sinusoidal scale, we divided the 
measured power spectra by a smooth 'wiggle-free' reference 
spectrum. For our purposes, this is the 'no- wiggles' spectrum 
of Eisenstein & Hu (1998) (see Blake & Glazebrook (2003); 
with real data, additional smooth polynomial terms can be 
fitted to remove any residual shape). We therefore do not 
utilize any information encoded by the shape of the power 
spectrum. The purpose of our philosophy is to maintain 
maximum independence from models and systematic effects: 
the shape of P(k) may be subject to smooth broad-band 
systematic tilts induced by such effects as complex biasing 
schemes, a running primordial spectral index, and redshift- 
space distortions. For the acoustic oscillations analysis, the 
power spectrum is measured in bins of width Ak — 0.01 h 
Mpc -1 . Plots of simulated power spectra divided by refer- 
ence spectra for different survey configurations are displayed 
in Figure 2. 

The resulting sinusoidal modulation for each realiza- 
tion is fitted with a simple 2-parameter empirical formula 
describing a decaying sinusoid, i.e. Equation 3 from Blake 
& Glazebrook 2003: 



P(k) 



= 1 + A k exp 



k 



0.1 h Mpc" 



/"2nk\ . , 

sin (id- (9) 



For each realization, we recorded (i) the best-fitting char- 
acteristic scale kA , (ii) the value of the chi-squared statistic 
for the 'no- wiggles' model (i.e. Equation 9 with A = 0), 
and (iii) the value of the chi-squared statistic for 



Xno — wi 

the best-fitting 'wiggles' model, xt 
was defined in the usual manner: 



The x statistic 



10 
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Figure 3. Histogram of values of P ro i (denned by Equation 11) 
for 400 Monte Carlo realizations of a survey with ri; m = 21.5 
and (to = 0.03, illustrating the skewed distribution of probabil- 
ities (note the logarithmic x-axis). The mean value of P re i (as 
plotted in Figure 4) is indicated by the vertical dashed line; real- 
izations possess less confident detections of acoustic oscillations 
than implied by this mean in just 19 per cent of cases. 



E 



Pohs{ki) -Pmodel(^i) 

a P (ki) 



(10) 



Our flat-sky approximation implies that the off-diagonal 
terms of the covariance matrix (i.e. correlations between ad- 
jacent power spectrum bins) are consistent with zero; this 
was explicitly tested by computing full covariance matrices 
for a test case. 

A simple relative probability of the 'no-wiggles' model 
and 'wiggles' model can be defined by 



P rc i = exp [-(xL- 



Xwig- 



0/2] 



(11) 



(but see below for further discussion). We note that the dis- 
tribution of values of P re i across Monte Carlo realizations of 
the Universe is far from symmetric, as illustrated by Figure 
3. In a Universe falling at the 50th percentile of the ensemble, 
the relative probability of the 'no-wiggles' model compared 
to the 'wiggles' model would be significantly lower than the 
mean of the distribution. 

As our initial assessment of the confidence of detection 
of the acoustic oscillations we considered the average value 
of the quantity P rc i defined by Equation 11 over the Monte 
Carlo realizations. We converted this into a probability for 
the 'no- wiggles' model by using P re i = P no -wig/-Pwig and 
fno-wig + fwig = 1. Contours of Pno-wig (expressed as a 
rejection 'number of sigmas' for a Gaussian distribution) are 
displayed in Figure 4 in the parameter space of (rn m , oo). In 
order to obtain a 3-cr detection confidence of 99.7 per cent 
(fno-wig = 3 x 10 ) we require a survey with parameters 
such that 



o-o < (n im - 19.5) x 0.01. 



(12) 



As an alternative method of quantifying the 'confidence 
of detection' of acoustic oscillations (i.e. the probability that 
A in Equation 9) we considered the following Bayesian 
approach. We placed a uniform prior Prior(A) on the value 
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Figure 4. Contours of the detection confidence of acoustic oscil- 
lations defined by the average value of P no — wig = Prcl/(1 + Prel) 
(Equation fl) for photometric-redshift imaging surveys with 
varying magnitude threshold ry lIa and photometric-redshift er- 
ror parameter crn. The probabilities are expressed as a rejection 
'number of sigmas' for a Gaussian distribution. We use a survey 
area of 10,000 deg 2 . 




1 2 3 

Acoustic oscillations amplitude A 

Figure 5. Probability distributions (as defined by Equation 14) 
of amplitude A for the first 5 Monte Carlo realizations of a sur- 
vey with parameters ri; m = 21.5 and <rn = 0.03. The curves are 
normalized such that P(Ab os t) = 1. The range of intercepts at 
A = illustrates the distribution of values for P no — wig across the 
Monte Carlo realizations. 



of A: 
Prior(A) 



A w id 

= elsewhere 



Amin < A < A mm + A 



(13) 



We chose A m i n = and A wid = 3 « 2Aacdm, to be 
conservative. We assumed prior knowledge of the acoustic 
wavescale Ica = 2-k/s, where s is the value of the sound 
horizon at recombination, known very accurately from lin- 
ear CMB physics (e.g. Eisenstein & Hu 1998, Equation 26). 
For an individual power spectrum realization, the probabil- 
ity density as a function of amplitude A is 



P(A)cxexp(-x72) 



(14) 



where \ 2 is the value of the chi-squared statistic of the fit of 
Equation 9 to the data of that realization (with kA = 2n/s). 
Figure 5 displays curves of P(A) against A for the first few 
Monte Carlo realizations of a simulated survey with rn m — 
21.5 and cr = 0.03. 

According to Bayesian statistics, the relative probabil- 
ity of a 'no-wiggles' and 'wiggles' model for one realization 
is: 

Pno — wig _ P(A — 0) (15) 



(r 



P(A) Prior(A) dA 



In the numerator of this expression, a 5-function prior cen- 
tred at A — has been integrated over. Figure 6 plots the 
average value over the Monte Carlo realizations of the quan- 
tity Pno-wig defined by Equation 15 (expressed as a rejec- 
tion 'number of sigmas' for a Gaussian distribution) in the 
parameter space of (nim,o"o)- Note that less confident de- 
tections of the acoustic oscillations are implied by using 
this prior, requiring surveys approximately half a magni- 
tude deeper for a 3-cr detection. This is reasonable because 
the probability density of the 'non-detection' model with 
A — is being compared to the average probability density 




r magnitude threshold 

Figure 6. Contours of the detection confidence of acoustic os- 
cillations defined by the average value of the Bayesian quantity 
Pno— wig (Equation 15) for photometric-redshift imaging surveys 
with varying magnitude threshold ri; m and photometric-redshift 
error parameter <rn. The probabilities are expressed as a rejection 
'number of sigmas' for a Gaussian distribution. We use a survey 
area of 10,000 deg 2 . 



of models with A ^ 0, rather than only to the best-fitting 
'detection' model. This serves to illustrate the critical role of 
priors in quantifying the 'confidence of detection'. Note that 
if we had widened our prior yet further by allowing a range 
in possible model acoustic oscillation scales kA then an even 
higher-performance survey would be required to achieve the 
same detection confidence. An alternative prior on the am- 
plitude for the 'wiggles' model would have been to use a 
delta function Prior( J 4) = S(A ~ Aacdm), which would pro- 
duce very similar results to using P re i (of equation 11). 
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Figure 7. Contours of the accuracy of determination of the char- 
acteristic scale kA for photomctric-redshift imaging surveys with 
varying magnitude threshold rn m and photometric-redshift error 
parameter <tq- We use a survey area of 10,000 deg 2 . 



3.2 Accuracy of measurement of acoustic 
oscillations 

We can also use our simulations to quantify the accuracy 
with which the characteristic scale (i.e. standard ruler) can 
be recovered from a photometric redshift survey, as a func- 
tion of rn m and ao- This is easily obtained as the spread 
in best-fitting values of kA over the Monte Carlo realiza- 
tions. We defined this spread as half the difference between 
the 16th and 84th percentiles of the distribution of fitted 
wavelengths. This quantity is plotted in Figure 7 as a per- 
centage fractional precision AA^/Aa in the parameter space 
of ()"ii m , cto). The precision improves with both increasing 
rii m and decreasing ao, peaking at ~ 0.7 per cent for our 
highest-performance survey (rij m = 24, ao = 0.01). 

In cosmological terms, this precision is equal to the ac- 
curacy with which the quantity x(z e s)/s can be determined 
by the survey (where x is the co-ordinate distance to the ef- 
fective redshift of the survey and s is the value of the sound 
horizon at recombination). This may in turn be converted 
into confidence distributions for dark energy models (e.g. 
Seo & Eisenstein 2003; Amendola et al. 2004; Dolney et al. 
2004; Glazebrook & Blake 2005). 



3.3 Comparison with spectroscopic surveys 

It is of great interest to compare the confidence and accuracy 
of the acoustic oscillation measurement from putative pho- 
tometric redshift surveys with those resulting from future 
spectroscopic redshift surveys. We therefore created Monte 
Carlo power spectrum realizations of a grid of spectroscopic 
surveys, using the same techniques as our photometric sur- 
vey analysis. We varied the total survey area An (from 1000 
deg 2 to 10,000 deg 2 ) and the limiting magnitude threshold 
riim (from 18 to 24). 

For our spectroscopic survey analyses we assumed the 
same redshift distributions as a function of rn m listed in 
Table 1, although we note that a realistic spectroscopic sur- 
vey would more likely be directed at a sub-population such 




r magnitude threshold 

Figure 8. Contours of the detection confidence of acoustic oscil- 
lations defined by the average value of P no -wig = Prel/(1 + Prel) 
(Equation 11) for spectroscopic redshift surveys with varying 
magnitude threshold ri; m and survey area An- The probabilties 
are expressed as a rejection 'number of sigmas' for a Gaussian dis- 
tribution. This plot may be compared directly with Figure 4 for 
photometric redshift surveys (note the different ranges of x-axis). 



as star-forming galaxies with strong emission lines, which 
would be selected in a more complex manner than a simple 
magnitude cut. 

Our spectroscopic redshift power spectra were analyzed 
for acoustic oscillation measurement in an identical manner 
to the photometric redshift surveys. For purposes of com- 
parison we bin power spectra averaging over angles, and do 
not separate the results into tangential and radial compo- 
nents. Figure 8 displays the confidence of detection as a 
function of (An,rn m ), quantified by the value of Pno-wig 
in the same manner as Figure 4. A 3-a detection of the 
acoustic oscillations can be achieved by a spectroscopic sur- 
vey with parameters (An = 1000 deg 2 , riim = 22.5) or 
(An = 3000 deg 2 , r Um = 21). For comparison, an equivalent 
detection is yielded by a 10,000 deg 2 photometric redshift 
survey with parameters (ao = 0.05, m m = 23) or (ao = 0.01, 
Him = 20.5). Note that the confidences of detection listed 
here are more conservative (by up to a factor of four in 
terms of the number of standard deviations for a Gaussian 
distribution) than those which would result from a full fit of 
a ACDM model power spectrum, as discussed and compared 
in Section 5.1. 

Figure 9 displays the resulting accuracy of measurement 
of the characteristic acoustic scale; this plot may be com- 
pared directly with Figure 7. For example, a spectroscopic 
survey of depth rn m ss 22.5 over An ~ 1000 deg 2 will achieve 
a 2% measurement of the standard ruler (a similar precision 
is achieved by a 10,000 deg 2 photo-z survey with the same 
depth and redshift error ao = 0.03). 

Figures 10 and 11 continue the comparison of photomet- 
ric and spectroscopic surveys. In Figure 10 we plot the ratio 
of areas of photometric and spectroscopic surveys achiev- 
ing the same accuracy of standard ruler measurement for 
a fixed magnitude threshold common to both surveys. We 
assume an area of 10,000 square degrees for the photomet- 
ric redshift survey and vary the area of the spectroscopic 
redshift survey, although the results are expected to apply 
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r magnitude threshold 

Figure 9. Contours of the accuracy of determination of the char- 
acteristic scale kA for spectroscopic redshift surveys with varying 
magnitude threshold r\ im and survey area An- This plot may be 
compared directly with Figure 7 for photometric redshift surveys 
(note the different ranges of x-axis). 
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Figure 11. The magnitude threshold required by a 10,000 deg 2 
photometric redshift survey (with (To = 0.03) to match the stan- 
dard ruler accuracy of a grid of spectroscopic surveys with varying 
area and depth. 
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Figure 10. The area ratio of photometric and spectroscopic red- 
shift surveys achieving the same accuracy of standard ruler mea- 
surement for a fixed magnitude threshold. This factor is deter- 
mined by the photometric redshift precision. 



more generally. From Figure 10 we see that, for a photomet- 
ric redshift precision <7n = 0.03, the area ratio for a fixed 
magnitude threshold is about a factor of 12. This is simply 
understood by the requirement that the number of Fourier 
modes contributing to the power spectrum measurement 
(m oc k x ,maxAn oc An/cro) must be roughly equal in the two 
cases. For example, if <to = 0.03 then kx,max ~ 0.02 h Mpc - . 
However, for a spectroscopic survey k x , max = fcu n « 0.2 h 
Mpc -1 , thus the same number of modes m is delivered by 
a survey area An reduced by a factor of « 10. The rela- 
tion m oc An/ao ~ constant also explains the overall scaling 
An oc (Jo apparent in Figure 10. 

A comparison at common magnitude threshold is of 
course simplistic: for given observational resources, an imag- 
ing survey can readily probe to fainter magnitudes. There- 
fore, Figure 11 considers a grid of spectroscopic surveys 
(parameterized by An, Him) an d illustrates the magnitude 



depth required by a 10,000 deg 2 photometric survey with 
an = 0.03 to match the standard ruler accuracy. The con- 
tours (r p hoto = 21 — » 24) correspond to standard ruler ac- 
curacies in the range 5 — * 1.5 per cent (see Figure 7). 



4 MEASURING THE TURNOVER USING 
PHOTOMETRIC REDSHIFT SURVEYS 

According to standard cosmological theory the clustering 
power spectrum should exhibit a 'turnover' (i.e. a maxi- 
mum in power) at a characteristic scale fcturn ~ 0.015 h 
Mpc -1 . The turnover arises because the primordial power 
spectrum laid down by cosmological inflation - hypothe- 
sized to be a featureless power law P pr im(fc) oc fc nBcalar , 
where n sca i ar ~ 1 - is suppressed by a 'transfer function' 
T(k) owing to radiation pressure in the radiation-dominated 
epoch. The resulting linear-theory power spectrum is derived 
as P(Jfc) oc P pri m(fc)T(fc) 2 where T(k) « 1 for k < fcturn 
and T(k) decreases towards zero for k > fctum with an ap- 
proximate limiting dependence k~ 2 . The characteristic scale 
fctum is equivalent to the co-moving horizon scale at matter- 
radiation equality. This is sensitive to the quantity f2 m /i 2 , 
since the larger the physical density of matter (f2 m /i 2 ), the 
earlier matter-radiation equality occurs, and suppression of 
growth due to radiation oscillations below the Jeans length 
does not have time to reach larger scales. Therefore the scale 
of the turnover is smaller and fct urn is larger. (Note that since 
we measure redshift and not distance, the a;— axis of the 
power spectrum plot is in units of h Mpc -1 and so the posi- 
tion of the turnover on a plot of P(k) against k/h depends on 
fl m h). Moreover, structure modes with wavelengths larger 
than the turnover scale are relatively unaffected by physics 
subsequent to inflation and potentially constitute a probe of 
the inflationary epoch. 

Detection of the turnover in the galaxy clustering pat- 
tern constitutes an interesting test of the cosmological 
paradigm. Its absence may imply either a failure of the 
standard cosmological theory, or the discovery of new large- 
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Figure 12. Mean values of the power spectrum P(k) averaged over the Monte Carlo realizations for various surveys, together with the 
standard deviation in each bin. In this plot we do not divide by any smooth reference power spectrum. The left-hand panel illustrates 
the variation of the obtained P(k) as a function of r Jim for <tq = 0.03; the right-hand panel displays the dependence on ctq for r lim = 22. 
Results for different surveys are offset for clarity. The dashed line tracks the value of fcx,max for the surveys in question. For a bin at 
k < fci.max, no further improvement in power spectrum precision is possible via more accurate values of arj, rather a survey must go 
deeper, mapping more cosmic volume and increasing the density of states in fc-space. Power spectra are only plotted for scales k < 0.05 h 
Mpc -1 . 



scale galaxy biasing mechanisms (e.g. Durrer et al. 2003; 
i.e., the turnover is more susceptible to systematic distor- 
tions than the acoustic oscillations). Successful definition of 
the turnover requires a survey possessing an extremely large 
volume to reduce the effect of cosmic variance. In addition, 
the number density of the tracer galaxies must be sufficient 
to suppress the shot noise contribution to the power spec- 
trum error. For example, quasi-stellar objects can easily be 
detected to high redshift but possess an inadequate number 
density to permit an experiment limited by cosmic variance 
(Outram et al. 2003). As a result, no survey has cleanly mea- 
sured the turnover yet (in a model-independent manner) . 

For a turnover detection experiment, spectroscopic- 
redshift accuracy is not required. The relevant scales are suf- 
ficiently large that equivalent information may be recovered 
from a photometric redshift survey, if the main contribu- 
tion to the photometric redshift errors is statistical and not 
systematic. For example, for a survey with <to = 0.03 and 
z c ff = 0.5, all Fourier modes with k x < fc^max ~ 0.02 h 
Mpc" 1 > fctum survive the radial damping. 

We now consider the detectability of the turnover and 
implied accuracy of determination of the characteristic scale 
fctum for a series of simulations varying nim and <jo. As be- 
fore, our starting point is the ensemble of power spectrum 
realizations obtained as described in Section 2.2. For the 
turnover analysis, the power spectrum was measured in bins 
of width Afc = 0.005 h Mpc" 1 . Because the turnover in P(k) 



wc 
l 



occurs at an approximate scale of fct urn ~ 0.015 h Mpc -1 
only utilize power spectrum modes with k < 0.04 h Mpc 
(i.e. 8 bins). Plots of simulated power spectra for different 
survey configurations are displayed in Figure 12. 

As with the acoustic oscillations analysis, the signif- 
icance of detection of the turnover depends on our prior 
assumptions. We again take a conservative approach, fitting 
our realized power spectra with a simple empirical parabolic 
turnover model characterized by four parameters: 



P(k) = P 
= Po 



1 - a 



1 



/ fc-fcp y 
V k J 



(k < fco) 

(k > fc ) (16) 



The free parameters are the turnover scale fco, the maxi- 
mum of the power spectrum Po, and the amplitudes of the 
parabolic decrease of power on either side of the maximum, 
a and f3. In this sense, a detection of the turnover is de- 
termined by finding a best-fitting value for a significantly 
greater than zero, and is governed solely by power spectrum 
modes at scales larger than the turnover scale (i.e. fc < fco). 
The requirement that P(k) > restricts the fitted param- 
eters to lie in the ranges Po>0, a<l,/3<l and we 
implement these conditions as strong priors in our fitting 
process. 

By analogy with the acoustic oscillations analysis, we 
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Figure 13. Contours of turnover detection confidence defined by 
the average value of P no — turn — ^rel 

/(I + P re i) (Equation 17) 
for photometric-redshift imaging surveys with varying magnitude 
threshold rii m and photometric-redshift error parameter cro. The 
probabilities are expressed as a rejection 'number of sigmas' for 
a Gaussian distribution. We use a survey area of 10,000 deg 2 . 

defined the confidence of turnover detection for a given 
power spectrum realization by comparing the chi-squared 
statistic for the best-fitting turnover model (i.e. Equation 
16) with that for the best-fitting 'no turnover' model, which 
we defined as Equation 16 with a set equal to zero: 

P r cl = exp [-(Xno-turn ~ Xturn-bost ) /2] (17) 

Figure 13 plots the average value of P no -tum = frel/(l+-Prei) 
over the Monte Carlo realizations (expressed as a rejection 
'number of sigmas' for a Gaussian distribution) in the pa- 
rameter space of (rum, cro). I n our highest-performance sur- 
vey (ri im = 24, cro = 0.01), the turnover is detected with 
~ 99.5 per cent confidence (Pno-turn = 0.005). A 2-cr de- 
tection in the mean realization requires rii m ~ 22.5. We 
note that, provided cro < 0.04, the detection confidence is 
independent of the photometric redshift accuracy because 
all Fourier modes beyond the turnover are retained in our 
analysis. 

Figure 14 quantifies the accuracy with which the char- 
acteristic turnover scale fcturn can be recovered from a pho- 
tometric redshift survey, as a function of rii m and cro, using 
the same technique as for the characteristic acoustic oscilla- 
tion scale in Section 3. In the best case we considered, the 
turnover scale can be measured with a precision of « 12 per 
cent. This is considerably poorer than the measurement ac- 
curacy of the acoustic oscillations scale, owing to the broad- 
ness of the turnover and the vastly fewer Fourier modes 
available at the relevant scales. This observation could in 
principle yield a 12% measurement of f2 m /i. This in itself is 
not particularly competitive with other techniques, but per- 
forming the analysis in this model-independent way tests the 
fundamental assumptions made in the standard analyses, 
and isolates the possible influences of relevant systematic 
effects. For example, it may indicate that scale-dependent 
biasing occurs on large scales (e.g. Dekel & Rees 1987) which 
would constitute a critical observation in the field of galaxy 
formation. Furthermore, precise measurements of the large- 




r magnitude threshold 

Figure 14. Contours of the accuracy of determination of 
the characteristic turnover scale fctum for photometric-redshift 
imaging surveys with varying magnitude threshold rn m and 
photometric-redshift error parameter cro. We use a survey area 
of 10,000 deg 2 . 

scale clustering pattern may unveil currently-unknown sig- 
natures of inflation or of non-Gaussianity (e.g. Martin & 
Ringeval 2004) pointing to a new cosmological paradigm. 

Unlike the case of acoustic oscillations measurement, 
we do not make a comparison with spectroscopic redshift 
surveys for our turnover detection experiment, because all 
relevant Fourier modes are retained by a photometric red- 
shift survey and therefore the results would be unchanged 
if perfect redshifts were known. Thus photometric redshift 
surveys will always be more efficient for turnover detection 
if the systematic errors can be controlled. 



5 MEASUREMENTS OF THE 

COSMOLOGICAL PARAMETERS 

In this Section we adopt a more model-dependent analy- 
sis approach to that performed in Sections 3 and 4, assum- 
ing the full theoretical framework of ACDM transfer func- 
tions and calculating how our simulated measurements of 
the galaxy power spectrum from photometric redshift sur- 
veys can be used to constrain cosmological parameters more 
tightly. Our investigation here is thus independent, but com- 
plementary, to the results presented earlier, and indicates 
how the cosmological conclusions are tightened by the in- 
corporation of more model assumptions. 

We assume linear biasing, i.e. that the galaxy power 
spectrum is a constant multiple of the matter power spec- 
trum, and marginalize over this parameter with a flat prior. 
We assume that the bias parameter does not evolve with 
redshift, and discuss the effect of this approximation in Sec- 
tion 6. 

In order to search parameter space we use 
Markov-Chain Monte Carlo (MCMC) sampling using the 
Metropolis-Hastings algorithm. In short, an MCMC 'chain' 
is made up of a list of 'samples' (coordinates in parame- 
ter space) which are obtained from performing trial likeli- 
hood evaluations Pr(x). A new sample at position Xi+i is 
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accepted with a probability min(Pr(xi+i)/Pr(xi), 1). The 
difficulty lies in suggesting good trial positions; we use the 
latest version of CosmoMC (see Lewis & Bridle 2002 and 
http://cosmologist.info/cosmomc for more information) 
which uses CAMB (Lewis, Challinor & Lasenby 2000) to 
calculate CMB and matter power spectra. We ran MCMC 
chains for the least powerful photometric redshift survey and 
used importance sampling to find parameter constraints for 
the better surveys. 

We used flat priors on the CMB parameters fib ft 2 , fi c /i 2 , 
dpcuk, t, n B , n r un and log 10(A S ). # pca k is used instead of the 
Hubble constant h because it renders the MCMC method 
more efficient when CMB data is included; it is defined by 
approximate formulae which given the CMB first peak po- 
sition in terms of the other cosmological parameters. There- 
fore for a given set of cosmological parameters, # pca k can be 
converted into h, and vice versa. For the simulations in this 
Section we take as our fiducial parameters those from the 
abstract of Spergel et al. (2003): h = 0.72, fi b /i 2 = 0.024, 
fi m /i 2 = 0.14, r = 0.16, n s = 0.99. The widths of the priors 
are chosen to be sufficiently large that they have no influ- 
ence on the results. We assume adiabatic initial conditions 
with a negligible tensor contribution. 

5.1 Measurements of the baryon fraction 

First we considered constraints on the baryon fraction re- 
sulting from galaxy surveys alone, as a function of Q m h. We 
fixed all other cosmological parameters at their input values 
for ease of comparison with current galaxy power spectrum 
analyses (e.g. Cole et al. 2005) . We derived results for simu- 
lated power spectra from the SDSS main spectroscopic sur- 
vey (using the survey parameters listed in Section 2.4) and 
for an SDSS photometric redshift imaging survey (assuming 
(To = 0.03, nim = 21) using sky area 10,000 deg 2 in both 
cases. These contours are plotted as the darker lines in Fig- 
ure 15. Photometric redshifts from the SDSS imaging survey 
would produce tighter parameter constraints than the final 
SDSS spectroscopic survey, if redshift accuracy ero = 0.03 
to a magnitude limit of nim = 21 could be achieved. This 
is a challenging requirement, but may be approachable by 
selecting Luminous Red Galaxies. 

In order to quantify the contribution of the overall shape 
to the detection of fib/fi m we added in as free parameters 
the Hubble constant and the primordial power spectrum tilt 
parameter n$. We applied an HSTKP prior on the Hub- 
ble constant, otherwise the Hubble constant had significant 
probability above 100 km s -1 Mpc -1 . The new constraints 
from the projected SDSS spectroscopic and photometric red- 
shift surveys are shown by the lighter contours in Figure 15. 
This is a more rigorous test of detection of acoustic oscil- 
lations - hence the slightly wider contours. The effect of 
relaxing the assumptions on n s and Ho is small when the 
constraints are weak, since the prior that fib > already 
limits the maximum error bar. 

We quantify the constraint on the baryon density by 
calculating the probability as a function of fib /fim marginal- 
ized over fi m /i and any other parameters. We then find the 
error bar by halving the distance between equi-probable lim- 
its containing 68 per cent of the probability and quote a 
'number of sigma' by dividing the fiducial value of fib/fi m 
by this error bar. The 'number of sigma' for our full SDSS 
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Figure 15. Contours in the parameter space of (£7b/^m, f2 m /i) 
derived from galaxy power spectrum measurements alone. To gen- 
erate the red and black (darker) lines we fixed other cosmological 
parameters at their input values. The dashed red (darker) contour 
is for our projected SDSS spectroscopic survey and the inner solid 
black (darker) contour is for a photometric survey with gq = 0.03 
and r < 21. For the cyan and green (lighter) lines Ho and n B are 
marginalized over using an HST Key Project prior of 72 ± 8 km 
s — 1 Mpc -1 . Only 68% confidence contours are shown, for clarity. 



spectroscopic survey simulation decreases from 3.5 to 2.3 on 
allowing Ho to vary within the HSTKP prior and completely 
freeing the spectral index n s . For the example photometric 
redshift survey shown here, it changes from 9 to 5. 

We note that the constraints on the baryon density from 
our analysis are rather weak from the main SDSS survey, 
despite the fact that it is larger than the full 2dFGRS sur- 
vey for which the baryon density has already been detected 
with better confidence. This is because we have assumed a 
relatively conservative value for the maximum wavenumber 
fitted by the linear power spectrum, fc max = 0.11 h Mpc -1 
(see Section 2.2). The exact results are quite sensitive to 
this value. In practice, experimental teams may choose to 
use a larger value of fc max , obtaining tighter contours than 
those displayed in Figure 15, but increasing their sensitiv- 
ity to the systematic uncertainties of modelling the quasi- 
linear regime. For comparison with existing measurements 
(e.g. Cole et al. 2005) and for maximum contrast with the 
'model independent' sections, we fix Ho and n s for the re- 
mainder of this Section. 

Figure 16 displays how the detection confidence of 
fib/fim depends on the parameters of a general set of future 
photometric redshift surveys. We conclude that the signif- 
icance of measurement of a non-zero value for fib/fi m will 
shortly be greatly improved by the use of photometric red- 
shifts. For a detection with 6-a confidence we require a sur- 
vey with rii m m 20 and ao < 0.04. 

Equivalent confidences are shown in Figure 17 for a 
general set of future deeper spectroscopic redshift surveys. 
Clearly for a given magnitude limit, the area required to 
achieve a detection of fib/fi m 7^ is smaller than the 10,000 
deg 2 used for the photometric survey simulations with the 
same magnitude limit, however to survey this area is signif- 
icantly more costly. A quantitative comparison between the 
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r magnitude threshold 

Figure 16. Input value of f2b/n m divided by the 68 per cent 
confidence error margin for O^/^m (loosely speaking, the 'num- 
ber of a of detection') for a range of photometric redshift surveys 
with varying magnitude threshold r and photometric-redshift er- 
ror parameter oq. The shape parameter fl m h is varied, as for the 
darker contours in Figure 15. We use a survey area of 10,000 deg 2 . 
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r magnitude threshold 

Figure 17. Input value of fl^/Qm divided by the 68 per cent con- 
fidence error margin for £7|->/f2 m (loosely speaking, the 'number of 
a of detection') for a range of spectroscopic redshift surveys with 
varying magnitude threshold r and area. The shape parameter 
Q m /i is varied, as for the darker contours in Figure 15. 

results for spectroscopic and photometric redshift surveys is 
shown in Figure 18. This is derived from the previous two fig- 
ures by calculating the area of spectroscopic survey required 
to obtain the same detection confidence as the 10,000 square 
degree imaging survey, for each magnitude limit and photo- 
metric redshift error. From this figure it can be seen that 
the factor in area required to make up for the photometric 
redshift uncertainties is about 12(cto/0.03) and roughly in- 
dependent of magnitude limit, in good agreement with the 
'model independent' analysis. 

The value of fib/fim is connected to the presence of 
the acoustic oscillations in the matter power spectrum, and 
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Figure 18. Area factor required to obtain the same accuracy of 
measurement of Q^/Qm for a photometric redshift survey as for 
a spectroscopic survey. 



thus, as discussed in Section 3, a high-significance measure- 
ment of a non-zero value for Qb/Q m could be considered 
as a potential 'detection' of the existence of these acoustic 
features. Our treatment in the current Section is thus more 
model-dependent than and complementary to our analysis 
method of Section 3. In the present Section we implicitly 
assume as a prior the entire ACDM framework in which the 
acoustic oscillation position, amplitude and matter power 
spectrum shape are intrinsically linked; whereas in Section 3 
we adopted a more conservative approach, simply fitting a 
modified sinusoidal function to the simulated data. 

Therefore it is not surprising that our detection con- 
fidences for fib/^m 7^ are somewhat tighter than those 
for a 'model-independent' detection of acoustic oscillations. 
Roughly, the 8-rr lines in Figure 16 lie on the 2-a lines of Fig- 
ure 4. Generally the number of a is a factor of four larger 
for the 'model-dependent' fit with h and n s held fixed. As 
illustrated by Figure 15, if we instead marginalize over the 
Hubble constant, with an HSTKP prior and free n s , then 
the measurements are much less accurate. The r < 21, 
(Jo = 0.03 survey constraint weakens from 9-a to 5-ct on 
freeing h and n s in this way, compared to 2.2-<r for the fully 
'model-independent' fit. 

The two different analysis methods presented in this 
study (i.e. Sections 3 and 5) are analogous to those used to 
detect the CMB EE polarization signal (DASI, Kovac et al. 
2002; CBI, Readhead et al. 2004). For the first detection of 
a non-zero signal, the DASI team used a 'template' for the 
shape of the EE power spectrum, taken from the ACDM 
model that best fit the CMB TT data. By contrast, the CBI 
team were the first to detect the phase and amplitude of 
the EE polarization signal using a model-independent sine- 
wave fit, and found that the inferred values were consistent 
with a ACDM model. In order to test the framework of 
the ACDM model we argue that a precise measurement of 
baryon wiggles in the matter power spectrum using a model- 
independent fit will be an enormous break-through; this be- 
ing done, more assumptions can then be made to extract 
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the most accurate possible measurements of the cosmologi- 
cal parameters. 

5.2 Measurements of the running spectral index 

A vital role for measurements of the galaxy power spec- 
trum is to break the parameter degeneracies inherent in the 
Cosmic Microwave Background (CMB) anisotropies. Using 
a standard six-parameter ACDM model, the WMAP satel- 
lite measurements can be readily converted to constraints 
on the matter power spectrum (e.g. Tegmark 2003), which 
are sometimes erroneously interpreted to mean that there is 
no need for galaxy redshift surveys. However, the most im- 
portant current questions in cosmology, such as the quest to 
quantify the properties of the dark energy and of inflation, 
demand that the simplest cosmological model be extended 
to encompass additional parameters, such as a time- varying 
dark energy equation-of-state and a more general model for 
the primordial power spectrum of mass fluctuations. In these 
cases the degeneracies inherent in the CMB become insuper- 
able, and high-quality additional data is essential. 

In the standard six-parameter cosmological model, it 
is assumed that the scalar perturbations have a power-law 
power spectrum parameterized by a single spectral index n B . 
However, the simplest inflationary models predict that this 
index should exhibit a slight dependence on scale, often pa- 
rameterized n run , such that the primordial power spectrum 
assumes the form 

P{k) cx (fc/fc ) ns+nrun log (k/ko) (18) 

where according to standard inflationary models, n run ~ 
0.002. This parameter n run has been the subject of much 
recent debate due to the apparent detection of a non-zero 
value by the WMAP team (Spergel et al. 2003) at n Iun = 
-0.031 ± 0.024 from WMAP, 2dFGRS, ACBAR and CBI. 
Therefore we include 

TT-run £111 additional parameter in our 
analysis, expanding the total number of fitted parameters to 
seven (noting that the presence of a running spectral index is 
potentially degenerate with that of a scale-dependent bias). 

Figure 19 plots the resulting constraints on the primor- 
dial power spectrum parameters of Equation 18 - n s and 
n run - when the WMAP CMB data is combined with our 
simulated 10,000 deg 2 SDSS main spectroscopic and pho- 
tometric redshift (ao = 0.03) surveys. The more accurate 
power spectrum measurements yielded by the photomet- 
ric redshift survey helps break the degeneracy between the 
scalar spectral index n s and the running spectral index n run 
(although the exact direction of the degeneracy is deter- 
mined by the pivot scale in Equation 18, for which we use 
k = 0.05 Mpc" 1 ). 

Figure 20 displays the (1-c) accuracy of measurement 
of n run for a general set of future photometric redshift sur- 
veys. For comparison, our projection for the full SDSS spec- 
troscopic sample is ±0.037. The constraint would be slightly 
tighter if more CMB data were included in the analysis. The 
error in n Iun is halved for our best possible imaging survey 
case (rn m = 24, ao = 0.01), bringing the limit in between 
that predicted by slow-roll inflation and that indicated by 
the WMAP first-year results. 

In Figure 21 we show the equivalent constraints for a 
general set of future spectroscopic redshift surveys. We see 
that to achieve an error bar <5n run = 0.03, a spectroscopic 
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Figure 19. Confidence contours in the parameter space of (n run , 
n s ) resulting from our simulated SDSS spectroscopic survey 
(outer contours) and from a simulated SDSS photometric red- 
shift survey with parameters ru m = 21 and co = 0.03 (inner 
contours). In each case, the galaxy power spectrum data is com- 
bined with WMAP CMB data and we marginalized over the other 
five cosmological parameters. 




r magnitude threshold 

Figure 20. 68 per cent confidence error margin of n run for a 
range of photometric redshift surveys with varying magnitude 
threshold r and photometric-redshift error parameter ctq. In all 
cases we combined with WMAP data and marginalized over the 
remaining cosmological parameters Hq, fi b /i 2 , Q m h 2 , t, erg and 
n s . 

survey requires about 15 per cent of the area of a photomet- 
ric redshift survey with ao = 0.03 to the same magnitude 
limit. 

As discussed above, detection of the matter power spec- 
trum shape on the largest scales is important because it has 
been unchanged since inflation. In Figure 22 we indicate 
the range of model matter power spectra permitted by the 
WMAP data (Verde et al. 2003, Hinshaw et al. 2003, Kogut 
et al. 2003) plus our fiducial SDSS photometric redshift sur- 
vey, for a seven-parameter ACDM model (i.e. including a 



© 0000 RAS, MNRAS 000, 000-000 



Cosmology with photometric redshift surveys 17 




magnitude threshold 



Figure 21. 68 per cent confidence error margin of n run for a 
range of spectroscopic redshift surveys with varying magnitude 
threshold r and area. In all cases we combined with WMAP data 
and marginalized over the remaining cosmological parameters Ho , 
Vt^h 2 , f2 m /i 2 , t, <rg and n s , as in Figure 20. 




Figure 22. The range of model matter power spectra allowed by 
the WMAP data combined with a photometric redshift survey 
with parameters erg = 0.03 and r < 21, colour-coded according 
to the value of n Tun . The plotted data points display the simu- 
lated photometric redshift power spectrum. Seven cosmological 
parameters are allowed to vary: Ho, SJt/i 2 , fi m /i 2 , t, as, n a and 
n run . Since we marginalize over the linear bias parameter then for 
display purposes the power spectra are normalized to go through 
the highest k data point. 

free parameter n Iun ). This shows that despite the improve- 
ments in power spectrum precision, there is still some free- 
dom in our knowledge of the matter power spectrum on large 
scales. 

Further relaxations of the post-inflation assumptions 
will heighten the importance of matter power spectrum in- 
formation. Whilst the increasing amount of CMB polariza- 
tion information will help to improve constraints, the range 



of possible models may be widened even further. For exam- 
ple: in the above we have assumed that the perturbations 
are adiabatic, with a negligible tensor contribution. More- 
over, in addition to adding tensors, a number of isocurvature 
modes are possible, along with freedom in their spatial cor- 
relations (e.g. Bucher et al. 2004) . In addition each mode has 
a power spectrum that could be relaxed from the power law 
and gentle running forms assumed here. We would ideally 
like to reconstruct these power spectra, or equivalently the 
inflationary potentials. For example a phase transition dur- 
ing inflation can cause a step- like feature in the scalar power 
spectrum, and trans-Planckian effects can induce 'ringing'. 
This additional freedom can be constrained effectively by 
combining CMB and large-scale structure information (e.g. 
Mukherjee & Wang 2003; Bridle et al. 2003). 



6 ASSESSING OUR APPROXIMATIONS 

We now assess the effect of the most significant approxi- 
mations contained in our methodology for simulating the 
accuracy of power spectrum measurements (Section 2). 

6.1 Photometric redshift error distribution 

Our fiducial set of simulations assumed that the statisti- 
cal distribution of photometric redshift errors could be de- 
scribed by a Gaussian function characterized by a standard 
deviation (Equation 2). This spread can always be measured 
by obtaining spectra of a complete sub-sample of the imaged 
galaxies. 

However, a real flux-limited survey will inevitably con- 
tain a combination of different classes of galaxy with differ- 
ent intrinsic photometric-redshift scatters. For example, Lu- 
minous Red Galaxies (LRGs) have especially strong spectral 
breaks which yield improved photometric redshift precision 
compared to an average galaxy possessing the same redshift 
and r-band magnitude. 

We assessed the effect of such combinations via sev- 
eral test cases in which the photometric redshift error dis- 
tribution was modelled by a sum of two Gaussian functions 
with different widths (denoted <7i, 02) and relative ampli- 
tudes (denoted bi , 62 such that 61 + 62 = 1). Specifically, 
we assumed o\ = 0.01, cr 2 = (0.03,0.1) and b 2 = (0.1,0.25), 
such that we are investigating the effect of a minority sub- 
population with a significantly broader error distribution 
than the majority of galaxies. 

Our motivation is to verify that the single Gaussian 
error function in Section 2 does not yield over-optimistic 
results relative to a more realistic double-Gaussian model 
with the same standard deviation, i.e. to check that our 
analysis is conservative. Thus in each case we compared 
the fractional power spectrum precision resulting from the 
double Gaussian model with that of a single Gaussian er- 
ror distri bution with the same overall standard deviation 
(oca = \Jb\o\ + b 2 a 2 ). 

For the double Gaussian model we derived the value of 
fc^.max in Equation 5 by taking the value of a x correspond- 
ing to the tighter of the two Gaussians (i.e. the dominant 
galaxy population) . We determined by experiment that this 
was optimal compared to other possible choices, such as a x 
corresponding to the overall standard deviation a e g. 
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We found that the precision of power spectrum mea- 
surement was never degraded by the assumption of a double- 
Gaussian model, and in the case (<72 = 0.1, 62 = 0.1) was 
significantly improved compared to our fiducial predictions 
for a single-Gaussian model owing to the tighter core of the 
photometric redshift error distribution. In other words our 
previous analysis is indeed conservative: the presence in real 
survey data of a minority sub-population of galaxies with 
significantly poorer photometric redshift precision than the 
overall standard deviation will not result in poorer power 
spectrum measurements than those inferred by the equiva- 
lent single Gaussian function. 

We also experimented with the addition of an overall 
systematic offset z a s to the mean photometric redshift (rel- 
ative to the spectroscopic value) of the second of the two 
Gaussian functions. For the case (02 = 0.03, 62 = 0.1) the 
power spectrum precision was degraded by w 10% when 
z a = 0.05. We note that such systematic effects can always 
be identified by spectroscopic follow-up of a large enough 
complete sub-sample of the imaging survey. 

6.2 Bias model 

Inference of the clustering pattern from galaxy surveys is 
always subject to uncertainties associated with the bias 
model, i.e. the precise manner in which galaxy light traces 
the underlying mass fluctuations. In our initial simulations 
we assumed that this biasing scheme was simply linear; 
but in general the bias mechanism will be non- linear, scale- 
dependent, non-local and evolving with redshift. 

Observationally, linear scale-independent bias appears 
a good approximation for a wide range of galaxy types 
on the large scales discussed here (e.g. Peacock & Dodds 
1994): structure formation is still in the linear regime and 
the physics of individual galaxy formation should not be 
relevant. We note however that halo-dependent effects may 
become important for the most massive galaxies (Peacock 
& Smith 2000; Seljak 2000). 

Even if bias is scale-dependent on large scales, it would 
be very surprising if it induced oscillatory features in k- 
space liable to obscure the distinctive acoustic peaks and 
troughs. Our model-independent analysis of acoustic oscilla- 
tions detection in Section 3 should be robust to such system- 
atic broadband tilts in the galaxy power spectrum because 
the overall shape is divided out. Some authors have argued 
that scale-dependent bias on large scales should obscure the 
turnover (Durrer et al. 2003); if this property was confirmed 
observationally then it would undermine our turnover de- 
tection analysis of Section 4, but would tell us something 
interesting about galaxy formation. 

In addition, we assume that the biasing scheme does 
not depend on redshift. This is in conflict with observa- 
tional data: Lahav et al. (2002) assume a model in which 
the clustering of galaxies remains constant with time (in co- 
moving space) whilst the dark matter perturbations grow, 
implying b(z) = b /D(z) where D(z) is the linear growth 
factor. Alternatively, if galaxies are assumed to follow the 
cosmic flow and remain constant in number then an al- 
ternative functional form is produced (Fry 1996) in which 
the linear bias parameter evolves in proportion to its devi- 
ation from unity in the present-day universe according to 
b(z) = 1 + {b - l)/D(z). Our fiducial model is b(z) = 1, 



which is a special case of the Fry equation, and in fitting 
to the data we have assumed the rather simplistic form 
b(z) = 6 . 

We note that despite their differences, all of these mod- 
els have the same number of free parameters, namely one, 
60 ■ Therefore if we were to fit instead the Fry bias model 
to a simulation with b(z) = 1 then we would have the same 
amount of information available for extracting cosmological 
parameters. Clearly b(z) = 1 is not a solution of the 'con- 
stant galaxy clustering' (CGC) model of Lahav et al. (2002) , 
but relative to this model our simulation and detections are 
conservative, since if 60 ~ 1 today (as measured for the 
2dFGRS by Lahav et al. 2002 and Verde et al. 2002) and 
b(z) — b a /D(z) (CGC model) then the galaxy clustering at 
higher redshift would be stronger than in our simulation and 
thus the signal to noise available for acoustic oscillation and 
turnover detection would be increased. 

We also point out that: 

• For a power spectrum measurement limited by cosmic vari- 
ance (such as those discussed here), the fractional error in 
the power spectrum (e.g. as plotted in Figure 2) is inde- 
pendent of the overall amplitude of the power spectrum as 
determined by the bias model. 

• An incorrect assumption about the bias model would not 
shift the acoustic oscillations along the wavenumber axis; 
the features would remain in the same place and add up 
coherently when integrating over redshift. 

6.3 Curved-sky effects 

Our flat-sky approximation would not be a valid analysis 
method for a real 10,000 deg 2 imaging survey: we must either 
break such a survey into smaller sky areas, or utilize a spher- 
ical harmonics decomposition along the lines of Percival et 
al. (2004). However, given that the fractional power spec- 
trum accuracy of our simulated surveys is determined prin- 
cipally by the cosmic volume surveyed, our results should 
be representative. In order to verify this, we re-analyzed the 
SDSS main spectroscopic survey discussed in Section 2.4 by 
casting the volume in a flat-sky box rather than a conical 
geometry (this case is at very low redshift, hence is particu- 
larly sensitive to this flat-sky approximation). The resulting 
power spectrum precision agreed to within 10 per cent. 

6.4 Analysis in redshift slices 

Application of the flat-sky analysis method to a real photo- 
metric redshift survey (with the redshift errors assumed to 
smear the distribution along one axis only) must addition- 
ally involve the survey being split into redshift slices: any 
significant ( > 20%) variation of the angular diameter dis- 
tance across an analyzed box would cause a smoothing-out 
of the acoustic features. Again our results should be repre- 
sentative because, regardless of the number of redshift slices 
employed, the same number of independent Fourier modes 
are being utilized (splitting a survey box into N slices de- 
creases the fc-space density-of-states in each slice by a fac- 
tor of N, but this factor is recovered by averaging over the 
slices). In order to verify this point we re-analyzed a test 
case, dividing a broad redshift range into several flat-sky 
slices such that the volume contained within each slice was 
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equal to the original curved-sky volume (i.e. increasing the 
width of the flat-sky slice with redshift). We additionally 
varied the power spectum amplitude in each slice in accor- 
dance with the cosmological growth factor. We obtained the 
final survey power spectrum by averaging the power spec- 
tra measured for the individual slices with inverse- variance 
weighting. The resulting power spectrum precision (recov- 
ered from the Monte Carlo realizations as usual) agreed with 
our original single-slice analysis to within 5 per cent. 

6.5 Angular selection functions 

In the present study we have not considered any effects 
due to complex angular selection functions, assuming that 
our survey area is simply a uniform box on the sky. For a 
more complex survey geometry (e.g. SDSS DR3), the ob- 
served galaxy power spectrum is smeared out according to 
the Fourier transform of the observed sky area (and correla- 
tions between adjacent power spectrum modes are induced). 
This convolution will smooth out features such as acoustic 
oscillations, reducing their detectability, thus surveys with 
reasonably contiguous geometries are required. For example, 
the angular selection function of the 2dFGRS (Percival et al. 
2001) causes a smearing of the galaxy power spectrum which 
is most significant on large scales (i.e. small k) and almost 
a delta-function at k ~ 0.1 h Mpc -1 (Elgaroy, Gramman & 
Lahav 2000). Such a selection function would therefore not 
seriously degrade measurement of the baryon oscillations; 
detection of the turnover would be largely unaffected if the 
width of the smearing was smaller than that of the turnover. 

6.6 Dust extinction 

Extinction by Galactic dust could affect the completeness of 
a survey with a given magnitude limit, and also compromise 
the accuracy of the photometric redshifts. This first prob- 
lem has to be addressed for both spectroscopic and imaging 
surveys alike, although the greatly increased depth and ar- 
eas covered by the imaging surveys imply that this will be 
a larger problem relative to the small random errors. This 
issue was discussed in detail by Efstathiou & Moody (2001), 
and was found to be unimportant for their results but would 
in general affect the galaxy power spectrum on the largest 
scales. For the purposes of this paper we have assumed that 
these effects can be overcome. The second problem affects 
only the photometric redshift surveys, and addressing this 
potential problem satisfactorily is beyond the scope of this 
paper. However, we note that in principle it could be over- 
come if the photometric redshift algorithm could be cali- 
brated as a function of dust column density, for example 
using a number of narrow spectroscopic surveys ('training 
sets') spanning a range of Galactic extinction optical depths. 



7 PROSPECTS FOR ONGOING AND FUTURE 
PHOTOMETRIC REDSHIFT SURVEYS 

We now discuss our results in the context of existing or pro- 
posed imaging surveys that may be utilized for photometric 
redshift studies. Table 2 provides a list of such surveys, cov- 
ering both optical and near infra-red wavebands. 



7.1 Photometric redshift estimation techniques 

A variety of techniques have been proposed in the liter- 
ature for derivation of photometric redshifts from multi- 
colour photometry. The photometric redshift performance 
depends on the method used, together with a complex com- 
bination of the galaxy magnitude and redshift, the filter set, 
the signal-to-noise ratio of the photometry, and the type of 
galaxy spectrum (in general, redder objects yield more ac- 
curate photometric redshifts). A detailed simulation of this 
myriad of factors is beyond the scope of this paper; how- 
ever, some general conclusions may be inferred using the 
photometric redshift accuracies discussed in the literature. 

The simplest photometric redshift techniques employ a 
limited set of 'template' spectra corresponding to local ellip- 
tical, spiral and starburst galaxies (e.g. Hyper-Z; Bolzonella, 
Miralles & Pello 2000). These templates may be redshifted 
and fitted to observed galaxy colours, deriving a likelihood 
distribution for the galaxy redshift. This approach has been 
successfully used in many cases (e.g. the Hubble Deep Field 
North; Fernandez-Soto, Lanzetta & Yahil 1999) and subject 
to various modifications such as the incorporation of mag- 
nitude priors in a Bayesian framework (Benitez 2000) and 
the iterative improvement of the initial templates (Budavari 
et al. 2000). Disadvantages of the method include the need 
for spectro-photomctric calibration over a wide wavelength 
range, the difficulty of incorporating galaxy evolution with 
redshift, and the failure of the template set to encompass all 
possible classes of observed galaxy. 

The availability of a training set - spectroscopic red- 
shifts for a complete sub-sample of imaged galaxies - is very 
helpful for 'tuning' the photometric-redshift technique. In 
Table 3 we list a number of ongoing redshift surveys that 
could potentially be used for this purpose. In one possible 
application of the training set, the galaxy redshift is ex- 
pressed as a polynomial in the colours, and the coefficients 
are fitted using the training set (Connolly et al. 1995). Fur- 
ther improvements are possible if a different polynomial is 
adopted in each of a series of cells in colour-space (Csabai 
et al. 2003) . Alternatively, an artificial neural network may 
be trained to deliver similar information (Firth, Lahav & 
Somerville 2003; Collister & Lahav 2004). 

7.2 Existing surveys: SDSS 

The largest ongoing galaxy survey (Table 2) is the SDSS. 
Very recently, the spectroscopic component mapping Lu- 
minous Red Galaxies has been utilized to obtain the first 
convincing detection of the acoustic scale (Eisenstein et 
al. 2005). The inferred accuracy of the standard ruler in 
this study was £s 4% and the rejection significance of 
fib/^m = was 3.4-cr. These results agree well with 
our own simulation (see Section 2.4; we assume an ex- 
panded area of 10,000 deg 2 ) in which the average value of 
P rc i = exp[-(Xno-wi g - Xwig-best)/ 2 ] (Equation 11) over 
the Monte Carlo realizations is 0.035 (i.e. 2.1-cr for the 
'model-independent' method - although 46% of realizations 
perform better than 3-cr - with a full ACDM fit expected 
to improve these figures as noted in Section 5.1). Our simu- 
lated standard ruler precision is 2.5% (the difference being 
explained by a scaling with survey area of V 'An). 

Turning now to analyses of the SDSS imaging compo- 
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Table 2. Existing and proposed photometric redshift imaging surveys; we only list projects mapping « 1000 deg 2 or more. For all future 
surveys the survey parameters are largely illustrative. 



Survey 


Waveband 


Depth 


Area (deg 2 ) 


Start date 


SDSS DR3 


ugriz 


r < 20.8 


5282 


Released 


UKIDSS 


JHK 


K < 18.5 


7500 


2005 


KIDS 


u'g'r'iV 


r' < 24.2 


1700 


2005 


VISTA Wide 


YJHK 


K < 19.5, Y < 22.0, J < 21.2, H < 20.0 


3000 


2006 


VISTA Atlas 


JK 


K < 18.2, J < 20.2 


20000 


2006 


Pan-STARRS 


gRIZY 


R < 27.2 


1200 


2008 


CTIO DES 


griz 


r < 24.1, g < 24.6, i < 24.3, z < 23.9 


5000 


2009 


VISTA DarkCAM 


u'g'r'iV 


r <25 


~10000 


2009 


LSST 


rBgiz 


r < 26.5, B < 26.6, g < 26.5, i < 26, z < 25 


15000 


>2012 



Table 3. Existing and proposed spectroscopic redshift surveys. 



Survey 


Selection criteria 


Area (deg 2 ) 


No. galaxy redshifts 


Status 


CFRS 


I AB < 22.5 


~ 0.1 


591 


Released 


CNOC-2 


R < 21.5 


1.5 


6200 


Released 


0OMBO17 (Photozs only) 


R < 24 


0.25 


10,000 


Released 


SDSS DR3 


r < 17.1 


4188 


374767 


Released 


VVDS CDFS (Le Fevre et al. 2004) 


Iab < 24 


21x21.6 arcmin 2 


1599 


Released 


VVDS Deep 


Iab < 24 


1.3 


50,000 


Ongoing 


VVDS Wide 


I AB < 22.5 


16 


100,000 


Ongoing 


zCosmos 


Iab < 23 


2 


90,000 


Ongoing 


SDSS-2dF LRG 


i < 19.5 plus colours for 0.4 < z < 0.8 


300 


10,000 


Ongoing 


SDSS LRG 


i < 19.2 plus colours for 0.15 < z < 0.4 


5000 


75000 


Ongoing 


SDSS-II 


r < 17.1 


10000 


~ 10 6 


Start 2005 


DEEP2 


Rab < 24.1 plus colours for z > 0.7* 


3.5 


65,000 


Ongoing 


KAOS 


TBD 


~ 1000 


~ 10 6 


Proposed 


SKA 


TBD 


~ 30000 


~ 10 9 


Proposed 



nent, Csabai et al. (2003) applied a variety of photometric 
redshift techniques to the Early Data Release, determining 
an overall rms redshift scatter Sz « 0.03 for r < 18, rising 
to Sz « 0.1 at r « 21, by which point systematic redshift 
discrepancies due to large photometric errors have become 
important (i.e. the effective magnitude limit for reliable ap- 
plication of photometric redshifts may be rn m w 20). The 
VVDS survey databases would constitute a suitable spectro- 
scopic calibration set (see Table 3) for photometric redshift 
techniques that require it. As discussed in Section 6.1, the 
inevitable presence of a small fraction of interlopers with 
significantly larger redshift errors does not compromise the 
scientific results. 

Combining these results with our Figures 4 and 13 we 
conclude that a photometric redshift analysis of the en- 
tire SDSS imaging database may not succeed in detect- 
ing features in the galaxy power spectrum in the model- 
independent manner discussed in Sections 3 and 4. For ex- 
ample even if ao = 0.03 and r < 20.5 were possible then Fig- 
ures 4 and 13 would imply a wiggle 'detection' with 90 per 
cent confidence and a turnover 'detection' with 80 per cent 
confidence. However the prospects for using the ACDM fit 
in Figure 16 are better, implying a detection of n b /fi m 
with a significance of 8-cr (for fixed values of h and n s ). 

However, it is clear that certain sub-classes of galaxy 
perform significantly better regarding photometric redshifts: 



SDSS red galaxies yield an accuracy twice that of blue galax- 
ies (Csabai et al. 2003), and Luminous Red Galaxies (LRGs) 
permit a rms redshift precision of cro = 0.02 for z < 0.55 
(Padmanabhan et al. 2005). In addition, these LRGs in- 
habit massive dark matter haloes and are consequently bi- 
ased tracers of the large-scale structure (we assume a linear 
bias factor 6 = 2). The resulting amplification of the cluster- 
ing strength implies that a lower space density is required 
to yield a given power spectrum accuracy, although in some 
models of galaxy clustering the amplitude of acoustic oscil- 
lations is diluted for the most massive galaxies (e.g. Peacock 
& Smith 2000). 

The optimal imaging approach is therefore to analyze 
an LRG photometric redshift catalogue (a suitable spectro- 
scopic training set is the SDSS-2dF LRG redshift survey, see 
Padmanabhan et al. 2005). We simulated such a catalogue 
using a redshift interval 0.2 < z < 0.7, assuming a Gaussian 
redshift distribution peaking at z = 0.45 of standard devia- 
tion z = 0.1. We supposed LRGs could be selected from the 
photometric data with a surface density So = 100 deg -2 and 
photometric redshifts measured with an accuracy oq = 0.02 
(i.e. Sz — <7o(l + Zen) ~ 0.03; see Padmanabhan et al. 2005). 
In addition, LRGs are assigned a linear bias factor b = 2. 
Some of these assumptions may be optimistic, in particu- 
lar Padmanabhan et al. note that the photometric redshift 
accuracy degrades for z > 0.55. We quantified the confi- 
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dence of detection of acoustic oscillations for a 10,000 deg 2 
survey. The average value of P IC \ is 0.056, corresponding to 
a detection slightly less confident than but comparable to 
the LRG spectroscopic database, with 28% of realizations 
possessed a 'model-independent' confidence exceeding 3-cx. 
The simulated accuracy of the standard ruler measurement 
is 3.8%. For the turnover analysis, our simulated detection 
confidence is low: we derive an average value of P rc i = 0.33, 
with an accuracy for the turnover scale of 27%. Again, indi- 
vidual realizations may perform significantly better. 

We conclude that the SDSS imaging dataset has the 
potential to yield a marginal model-independent detection 
of acoustic oscillations using a sub-sample of red galaxies. 
However, we caution that the currently-available third SDSS 
data release (DR3), which covers a sky area w 5000 deg 2 , 
possesses a complicated angular window function which ren- 
ders this experiment more difficult: as discussed in Section 6, 
the observed power spectrum is a convolution of the underly- 
ing power spectrum and the survey window function; if this 
latter possesses a substantial width in Fourier space then the 
oscillatory signal will be smoothed and consequently harder 
to detect. 

7.3 Future surveys 

We now turn our attention to future surveys. We first note 
that the availability of near infra-red imaging to appropri- 
ate depths is extremely valuable for galaxies with redshifts 
z > 0.4. In this range, Bolzonella, Miralles & Pello (2000) 
derived a factor of 2 improvement in redshift accuracy when 
JHK photometry was added to the standard optical wave- 
bands. The combination of the future UKIDSS infra-red 
data (see Table 2) with the ongoing SDSS optical survey 
will therefore be very powerful, and should permit detection 
of the acoustic oscillations using the complete galaxy pop- 
ulation rather than special sub-classes. For example, if the 
photometric error were halved from ao — 0.03 to er = 0.015 
(assuming a magnitude threshold r = 20.5) then the model- 
independent wiggle detection limits in Figure 4 show a de- 
tection at 99 per cent confidence (as opposed to just over 90 
per cent without the infra-red data). 

It is clear from Figure 4 that in order to deliver a 
high-significance measurement of the acoustic peaks using 
photometric redshifts (together with a significant detection 
of the power spectrum turnover), a much deeper optical 
database is required (r ~ 24) over an area of several thou- 
sand deg 2 . Such a catalogue may first be provided by the 
CTIO Dark Energy Survey (DES; see Table 2). In this case, 
a photometric-redshift precision ao = 0.05 suffices for de- 
lineation of the acoustic oscillations and turnover, a red- 
shift accuracy which has been achieved in existing analy- 
ses of the Hubble Deep Field (e.g. Fernandez-Soto et al. 
1999). Furthermore, judged solely by the accuracy with 
which the power spectrum features can be mapped out, 
this photometric-redshift approach will be competitive when 
compared to the spectroscopic redshift surveys which will 
be contemporary to the DES (e.g. the KAOS proposal) - al- 
though we note that in terms of measuring the dark energy 
parameters, a spectroscopic survey yields critical additional 
information in the radial direction (namely, the Hubble con- 
stant at high redshift) which is forfeited by the photometric- 
redshift approach. 



The 'ultimate' photometric-redshift survey would cover 
an area An approximating the whole sky to a magnitude 
depth r ~ 26; surveys with the LSST will approach these 
specifications (see Table 2). If we assume that the redshift 
accuracy cannot exceed the limit er = 0.01, in this best 
case a competing spectroscopic survey would need to cover 
An/ 4: ~ 10000 deg 2 to produce comparable power spectrum 
constraints (i.e. trace the same number of Fourier structure 
modes). In optical wavebands, spectrographs with fields-of- 
view greatly exceeding 1 deg 2 are prohibitively expensive, 
but we note that redshift surveys for neutral hydrogen using 
a next-generation radio telescope such as the Square Kilo- 
metre Array would become competitive here (Blake et al. 
2004), provided that such a telescope was designed with an 
instantaneous field-of-view of order 100 deg 2 at frequency 
1.4 GHz. 

We note that ultimately the matter power spectrum de- 
rived from gravitational lensing will circumvent the assump- 
tions about bias necessary for the interpretation of galaxy 
surveys. Of course, lensing investigations can be performed 
with the identical imaging datasets discussed here, and in- 
deed the application of photometric redshifts will form an 
important part of that interpretation. However, due to the 
unknown intrinsic shapes of galaxies, many galaxies have to 
be averaged to obtain a cosmic shear signal. Therefore we 
will always be able to obtain constraints from the galaxy 
power spectrum that are vastly tighter than those inferred 
from the lensing power spectrum, albeit less reliable owing 
to the caveats regarding bias. Thus galaxy surveys are the 
best place to look for the first glimpses of any exciting new 
physics revealed by studies of large-scale structure. 



8 CONCLUSIONS 

We have used Monte Carlo techniques to estimate the mea- 
surement precision of the galaxy power spectrum achievable 
using photometric redshift imaging surveys with a variety of 
magnitude depths and photometric redshift accuracies. We 
have focussed on two main areas: 

• The 'model-independent' measurement of specific features 
in the galaxy power spectrum: the acoustic oscillations and 
the turnover. In particular, we have carefully quantified the 
statistical confidence with which these features may be de- 
tected, together with the accuracy with which the acoustic 
and turnover scales can be inferred. 

• The assumption of the ACDM paradigm to place tighter 
constraints on the baryon fraction f2b/fi m from galaxy sur- 
veys alone, and an evaluation of the additional constraints 
on the running of the spectral index of the primordial power 
spectrum as increasingly powerful galaxy surveys are com- 
bined with WMAP. 

We summarize our general findings as follows: 

• On linear-regime scales, a large photometric redshift sur- 
vey can provide competitive power spectrum measurements 
when compared to contemporaneous spectroscopic surveys. 
For example, given a magnitude threshold rn m = 23, a 
10,000 deg 2 photometric redshift survey with error parame- 
ter <To = 0.03 results in a confidence of detecting the acoustic 
oscillations that is comparable to a spectroscopic survey cov- 
ering 1000 deg 2 . (Albeit with the loss of radial information 
that is helpful for constraining dark energy models). 
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• We compare various different definitions of 'detection 
of acoustic osciliations', the most optimistic of which is a 
ACDM fit to the data (i.e. derivation of a baryon fraction 
significantly exceeding zero) which (for fixed values of h and 
n B ) yields detections with statistical significance approxi- 
mately a factor of four a's greater than our default 'model- 
independent' fit using no shape information. 

• Concerning the power spectrum turnover, the relevant 
scales are sufficiently large that photometric redshifts with 
precision better than ao — 0.04 retain information equiva- 
lent to spectroscopic surveys. An imaging depth r ~ 24 is 
required to deliver a 3-<r detection confidence. 

• Our results are robust against more complex photometric 
redshift error distributions such as double Gaussian models, 
if our parameter ao is interpreted as the overall rms error. 

In terms of realistic ongoing and future surveys: 

• The SDSS has already yielded a convincing detection of 
acoustic features in the spectroscopic survey of Luminous 
Red Galaxies. A marginal 'model- independent' detection 
may also follow from this sample. 

• Analysis of an LRG photometric redshift database selected 
directly from the imaging data may yield a measurement of 
comparable precision. 

• Assuming a ACDM model and calculating the number of 
standard deviations by which fib/fi m (marginalized over the 
shape parameter Q m h) exceeds zero shows that, for our sim- 
ulated full SDSS spectroscopic survey the baryons would be 
detected at 3.5-a, whereas for a photometric redshift sur- 
vey to magnitude depth r < 21 with Sz = 0.03(1 + z) the 
detection level is 9-er. We note that these analyses assume 
a conservative specification of the maximum wavenumber 
fcmax for which a linear power spectrum is fitted. 

• A high-significance analysis of the acoustic oscillations 
from photometric redshifts, and a model-independent de- 
tection of the power spectrum 'turnover', requires a signif- 
icantly deeper optical database (r ~ 24) over an area of 
several thousand deg 2 . 
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